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Abstract. This is the second paper of our investigation 
of the 0.5-2 keV soft X-ray luminosity function (SXLF) of 
active galactic nuclei (AGN) using results from ROSAT 
surveys of various depth. The large dynamic range of the 
combined sample, from shallow large-area ROSAT All-Sky 
Survey (RASS)-based samples to the satellite's deepest 
pointed observation on the Lockman Hole, enabled us to 
trace the behavior of the SXLF. While the first paper 



(Miyaji, Hasinger, Schmidt 2000 , hereafter paper I) em- 
phasized the global behavior of the SXLF, cosmological 
evolution and contribution to the soft X-ray background, 
this paper presents actual numerical values for practical 
use of our results. 

To present the binned SXLF, we have used an im- 
proved estimator, which is free from biases associated with 
the conventional ^ V^ -1 estimator. 
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1. Introduction 

The AGN/QSO luminosity function and its evolution with 
cosmic time are key observational quantities for under- 
standing the origin of and accretion history onto super- 
massive black holes, which are now believed to occupy the 
centers of most galaxies. Since X-ray emission is one of 
the prominent characters of the AGN activity, X-ray sur- 
veys are efficient means of sampling AGNs for luminosity 
function and evolution studies. An X-ray selected sample 
of AGNs is particularly useful because optical surveys of- 
ten use point-like morphology as a criterion for selecting 
AGNs (QSOs) among numerous other objects, and thus 
are likely to miss moderate-luminosity intermediate-high 



redshift AGNs embedded in their host galaxies. Also, ra- 
dio surveys sample only a minor population of AGNs. 

The Rontgen satellite (ROSAT) provided us with soft 
X-ray surveys with various depths, ranging from the 
ROSAT All-Sky Survey (RASS , Voges et al. |l999| ) to 
the ROSAT De ep Su rvey (RDS) on the Lockman Hole 
(Hasinger et al. 1998 ). Various optical identification pro- 
grams of the survey fields have been conducted and the 
combination of these now enabled us to construct the soft 
X-ray luminosity function (SXLF) as a function of red- 
shift. 

In paper I., we presented a number of global repre- 
sentations of the 0.5-2 keV SXLF and investigated the 
contribution to the soft X-ray background. We showed 
that our data are not consistent with the pure-luminosity 
evolution (PLE), contrary to the suggestions of a num- 
ber of previous analyses (e.g. Boyle et al. 1994; Jones 
et al. 1996). Instead, we find an excess of intcrmediate- 



redshift low-luminosity AGNs above the PLE case , som e 
sign of which was also recognized by Page et al. ( 1997 ). 
In view of this, we developed two versions of luminosity- 
dependent density evolution (LDDE1 and LDDE2) mod- 
els, which represent the observed data very well. An ex- 
trapolation of these two LDDE models below the faintest 
limit of the survey (S x = 1.7 x 10 -15 [ergs -1 cm -2 ]) yields 
significantly different predictions for fainter fluxes, brack- 
eting the range of a possible AGN contribution to the soft 
X-ray Background. Chandra (e.g. Mushotzky et al. [2000; 



Send offprint requests to: T. Miyaji 
* present address for TM 



Hornschemeier et al. gOOOJ) and XMM-Newton (gOOlJ) are 
probing much fainter sources and spectroscopic identifica- 
tions of these will eventually show which of these models is 
closer to the actual behavior of the AGN SXLF. However, 
because at least some of the faint X-ray sources are op- 
tically too faint for spectroscopic identification even with 
the largest telescopes, extending the XLF into such faint 
flux level may be difficult. 

In this second paper, we present practical and conve- 
nient expressions of the observed SXLF from the ROSAT 
surveys. We present our results mainly for the investi- 
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Table 1. R OS A T Surveys used in the Analysis 



Survey a 


crlim 
>->xl4 


Area 


No. of b 




KT^tergs^cm- 2 ] 


[deg 2 ] 


AGNs 


RBS 


« 250 


2.0 x 10 4 


216 


SA-N 


« 13 


685. 


130 


RIXOS 


3.0 


15. 


205 


NEP 


1.0 


0.21 


13 


UKD 


0.5 


0.16 


29 


RDS-Marano 


0.5 


0.20 


30 


RDS-LH 


0.17-0.9 


0.30 


68 



a Abbrevia tions - RBS: The R OSAT Bright Survey (Fis- 
cher et al. [1998| ; Schwope et al. |2000| ), SA-N: The S electe d 
Area-North (Zickgraf et al. |l997 ; Appenzeller et al. |l998j) , 



RIXOS: The R OSAT International X-ray Optical Survey (Ma- 
son et al. |l999|) , NEP: The North Ecliptic Pole Survey (Bower 
et al . |1996| ); UKD: The UK Deep Survey (McHardy et al. 
1998; ), RDS-Mara no: T he ROSAT Deep Survey - Marano field 
(Zamorani et al. |l999|, RDS-LH: The ROSAT Deep S urvey 



Lockman Hol e (H asinger et al 



Lehmann et al. 



1998; Schmidt et al. 1998 



199E). b Excluding AGNs with z < 0.015. 



gators who are interested in particular redshift regimes 
and/or comparing their models with observations. For this 
purpose, we show convenient analytical expressions in sev- 
eral redshift intervals separately. These give more accurate 
representations of the data in the redshift ranges of inter- 
est than those presented in paper I. We also tabulate the 
numerical values of the binned SXLF using an improved 
estimator. 

We use a Hubble constant Ho = 50 
[kms -1 Mpc~ 3 ]. The h 50 dependences are explicitly 
shown. We calculate the results with common sets of 
cosmological parameters: (f2 m ,f2A) =(1.0,0.0), (0.3,0.0). 
and (f2 m ,f2A) = (0.3,0.7). The symbol "Log" represents 
the base-10 logarithm. 

2. The summary of the sample 

We have used soft X-ray sources identified with AGNs 
with redshift information from a combination of ROSAT 
surveys in various depths/areas from a number of already 
published and unpublished sources. Detailed description 
of the definition of the sample, ROSAT countrate-to-flux 
conversion, and survey area are shown in paper I. The 
summary of the samples, which is a duplicate of Table 1. 
of paper I with updated references, is shown in Table |l|. 
The details of the nature and completeness of each sample 
were discussed in paper I. The limiting flux versus survey 
area relation were also shown in paper I. 

3. The SXLF estimation 

As in paper I, we present the SXLF in the observed 0.5-2 
keV band, i.e., in the 0.5(1 + z) — 2(1 + z) keV range in the 
object's rest frame. This is equivalent to assuming an en- 
ergy index of 1. Thus no K-correction was applied for our 



expressions presented here. The reasons for this choice are 
explained in detail in paper I. This choice is particularly 
important for this paper, which is intended to be used 
as observational constraints for population-synthesis-type 



models (e.g. Madau et al.1994; Comastri et al. 1995; Gilli 
et al. ( |1999[|200C| ); Miyaji et al. |2000D with various spectral 
assumptions. Because of that, it is more useful to provide 
quantities in a model-independent form rather than apply- 
ing a particular version of model-dependent K-corrections. 
By presenting in this manner, one can avoid the difficulty 
of reverse K-correcting and re-applying new K-corrections 
when the new results from Chandra and XMM provide 
better knowledge of the X-ray spectra of the population. 

3.1. Analytical expressions 

First, we find a smooth analytical function for each red- 
shift bin using a Maximum-likelihood fitting. The absolute 
goodness of the resulting expression can then be tested 
by one- and two-dimensional Kolgomorov-Smirnov tests 
(h ereaft er, 1D-KS and 2D-KS t ests r espectively; Press et 
alP 



1992 



1987]). See paper I for de- 
These fittings and 



Fasano & Franceschini 
tailed description of these methods, 
tests can be applied to unbinned data sets thus are free 
from artifacts and biases from binning. 

For an analytical expression, we use the smoothed two- 
power-law formula, as we did in paper I. Here, we fit the 
data in narrow redshift bins and thus evolution in each 
redshift bin is assumed to be a pure density evolution 
form: 



d$(L x ,z) 
d Log L x 



= A 



71 



72' 



1 + Z 

l + z c 



,(1) 



where z c is the central redshift of the bin. For the highest 
redshift bin where the "break" is not apparent, we have 
used a single power-law form by neglecting the first term 
in the square bracket in Eq.[l| 

The luminosity range of the fit is from Log L x =41.7 
to the maximum available luminosity in the sample. As 
shown below and in paper I, the SXLF below the minimum 
luminosity has a significant excess above the smooth ex- 
trapolation. This excess smoothly connects with the SXLF 



of the non-AGN population (e.g. Hasinger et al. 1999) and 
the X-ray emission may be significantly contaminated by 
non-AGN activities. 

The set of parameters which give the best fit for each 
redshift bin are shown in Table || along with the results 
of the ID- and 2D- KS tests (see the notes of the table). 
The parameter errors correspond to a likelihood change 
of 2.7 (90% confidence errors). In any case, Eq. [l] gives a 
statistically satisfactory expression for all redshift bins. 

3.2. TheJ^V- 1 Method 

The ^V^ 1 estimator, which is a gener alized version of 
the original V~^ x estimator (Schmidt 196§|) applied to 
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Table 2. Best- fit parameters for each redshift bin 



3 



2-range 


Zc 


N 






L% 


7i 


72 


P 


KS-prob b 












Ov) = 


(1.0,0.0) 








0.015-0.2 
0.2-0.4 
0.4-0.8 


0.1 
0.3 
0.6 


269 
113 

99 


(2.12 ±0.23) 
(3.72 ± 0.63) 
(1.55 ±0.28) 


■ 10~ b 

■ 10" 6 

■ 10~ 5 


0.60t{j;g 
0-891^° 
0.54irj 


u.oy_ 2 g 

u -°'-0.38 
ri qq+0.52 
u -°°-0.87 


o 1+0.4 
z - A -0.3 
9 r+0.4 
z '°-0.3 
9 9+O.3 
z - z -0.2 


4 00+ 2-53 
*- zz -2.61 
c 9C+3.48 
J - ZJ -3.51 
K qc + 2.29 


0.99,0.68,0.61 
0.93,0.56,0.72 
0.85,0.85,0.57 


0.8-1.6 


1.2 


135 


(1.75 ±0.27) 


■ 10" 5 




n 4n+ - 41 

u -^ u -0.53 


9 4+O.2 
z -*-0.2 


4.0711:12 


0.99,0.96,0.83 


1.6-2.3 


2.2 


44 


(3.75 ± 1.02) 


■ 10~ 5 


1.2(*) 


0.0(*) 


9 1 +0.2 
Z - A -0.1 


0(*) 


0.27,0.36,0.19 


2.3-4.6 


3.0 


25 


(3.93 ± 1.41) 


• 10 -5 


1.0(*) 




l q+0-2 


0(*) 


0.72,0.99,0.64 










(On, 


Ov) = 


(0.3,0.0) 








0.015-0.2 
0.2-0.4 
0.4-0.8 


0.1 
0.3 
0.6 


269 
113 
99 


(2.07 ±0.23) 
(3.33 ±0.56) 
(1.03 ±0.19) 


■ 10~ b 

■ 10~ 6 

■ 10~ 5 


0.591^ 

U.»O_ 49 

69+ 1 ' 07 


59+ u -^ 

u - Jy -0.30 

n c 7 +o.3i 

U - D ' -0.39 

n qy+0.50 

— 0.80 


9 1+0.4 
z - A -0.3 
9 4+O.4 
Z - 4 -0.3 

9 q+0.3 

z -°-0.2 


4-131^;^ 
5.3il 3 3 ;^ 

5.90l 2 ;^ 


0.99,0.60,0.54 
0.99,0.53,0.79 
0.97,0.81,0.57 


0.8-1.6 


1.2 


135 


(9.28 ± 1.44) 


■ 10~ 6 


z - A *-0.83 


42+ - 40 

u -^ z -0.52 


9 4+O.2 
z - 4 -0.2 


4-i3± 1 ;it 


0.98,0.97,0.72 


1.6-2.3 


2.2 


44 


(1.90 ±0.51) 


■ 10" 5 


1.8(*) 


0.0(*) 


9 1+0.2 
Z - A -0.1 


on 


0.19,0.41,0.15 


2.3-4.6 


3.0 


25 


(5.00 ± 1.80) 


■ 10" 5 


1.0(*) 




i q+0.2 
A - a -0.2 


on 


0.64,0.96,0.75 










(On, 


Qa) = 


(0.3,0.7) 








0.015-0.2 


0.1 


269 


(1.58 ±0.17) 


■ 10" b 


u - ' A -0.39 


62 +u -^ 

U - Oz -0.29 


9 1+0.4 

z ' A -0.3 


°- ' y -2.64 


0.99,0.54,0.61 


0.2-0.4 


0.3 


113 


(2.40 ± 0.41) 


■ 10~ 6 


1 nq+ 153 
A - uy -o.58 


fi7+ - 30 

U - D ' -0.39 


n 4+O.4 
Z - 4 -0.3 


i.95 


0.97,0.56,0.82 


0.4-0.8 
0.8-1.6 


0.6 
1.2 


99 
135 


(6.71 ± 1.22) 
(5.68 ±0.88) 


• 10" 6 
■ 10" 6 


2 69+ 2 ' 07 

Z - OM -1.04 


Sfi+ - 51 

43+0-39 


9 9+O.3 
z - z -0.2 

2 4 +0 - 2 

z '^-0.2 


5.69l 2 ; 28 7 

4-io±l;|! 


0.96,0.81,0.60 
0.98,0.97,0.74 


1.6-2.3 


2.2 


44 


(1.34 ±0.36) 


• 10 -5 


2.0(*) 


0.0(*) 


2 1+ 01 
z - A -o.i 


on 


0.14,0.40,0.16 


2.3-4.6 


3.0 


25 


(4.25 ± 1.53) 


• 10 -5 


1.0(*) 




i q+0-2 

A ' a -0.2 


on 


0.78,0.98,0.76 



Parameter values which have been fixed during the fit are labeled by '(*)'. a Units - A: [h$ Mpc -3 ], L„: [10 44 h^eig s -1 ]. b 
The three values are probabilities in two 1D-KS test for the distribution, L x , 1D-KS test for the z distribution and the 2D-KS 
test for the (L*,z) space respectively. 



a sample composed of subsamples of different depths (see 
paper I; Avni &; Bahcall 1980), has been widely used 



for binned luminosity functions (LF; we use the acronym 
LF when the discussion is not limited to the luminosity 
function in the X-ray band) in the literature. However, 
as discu ssed i n paper I (see also Wisotzki 1998 ; Page & 
Carrera 1999), using it for a binned LF estimator can 



cause significant biases, especially if the bin covers the flux 
range where the available solid angle of the survey changes 
rapidly as a function of flux. Also, the choice of the loca- 
tion in a Log L x bin with a non- negligible width at which 
the data point is plotted significantly changes the impres- 
sion of the plot. In Fig. 3 of paper I, however, we plotted 
the V& 1 estimates, because of the lack of a reasonable 
alternative at the time of writing that paper, with caveats 
on biases associated with the method. We note that the 
estimator can be used in an unbinned manner by consid- 
ering a set of delta- functions weighted by V~ x (or V~^) 
at the positions of samp le objects in the luminosity space 
(Schmidt & Green 1983 ), and this method is free from bi- 
ases mentioned above. While this unbinned method is a 
powerful tool to predict, e.g., the source counts, it does 
not provide practical means of plotting. In this paper, we 
have developed an improved estimator, which is explained 
in the next subsection. 



3.3. An improved estimator of the binned SXLF 

As an alternative to the ^V^ 1 method, we have devel- 
oped the following estimator for the binned LF, which is 
free fro m m ost biases unavoidable in the ^ V^ 1 method. 
In Sect. 



3.1 



we have found a smooth analytical function 
which describes the behavior of the SXLF in a given red- 
shift range. Having the best-fit smooth function, the esti- 
mated numerical value for the SXLF in a given bin in the 
(L x , z)-space is: 



d$ n 



d Log L x 



idl 



d Log L x 



(I* 



N: 



obs 



N 



mdl : 



(2) 



where L x ^ and Zi are the luminosity and redshift repre- 



sentative of the i-th bin, 



is the best-fit analytical 



d Log L x 

expression evaluated at this point, N° hs is the actual num- 
ber of AGNs observed in the i-th bin, and iV™ is the 
predicted number of AGNs in the bin from the best-fit 
analytical expression. Hereafter, we refer to Eq. || as the 
"JV obs /iV lndl estimat or" . N ote that the estimator proposed 
by Page & Carrera (1999) (hereafter PC) is a special case 



of Eq. H where d ^ L oc const (or 



d<S>" 



L x ). 



d L x k \ d Log L x 

Another advantage of this estimator over ~ V~ 1 is 
that exact errors at a given significance can be evaluated 
using Poisson statistics. One disadvantage of this estima- 
tor is that it is model-dependent, at least in principle. 
Since our analytical expressions are satisfactory represen- 
tations in any case and the estimator is not sensitive to 
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Simulation 



3. 



\ 
\ 

\ \ 



\ 



^: X! ^a 1 estimator 
O; N ob j/N rndl estimator 

n; Page & Carrera esLimaLor 



\ ^\ 
\ \ \ 
\ \ \ 

S \ N 
E\ © . 
\ \ 
\ \_ 

\ s 
\ 

\ 



42 43 44 45 

L °g k [ h 5n er e s l ] 

Fig. 1. The binned XLFs from a simulated sample us- 
ing three different estimators (symbols with error bars as 
labeled) are compared with the underlying "true" XLF 
represented by dashed lines. The vertical positions of the 
three different estimators have been shifted vertically for 
display. 



the details of the underlying model, the uncertainties due 
to the model dependence are practically negligible. 

In order to compare the goodness of the estimators, 
we performed simulations. Using the actual best-fit model 
for the 0.2 < z < 0.4 bin for (Q m , ft A ) = (1.0, 0.0), we 
generated a set of simulated AGNs. The number of sim- 
ulated AGNs are 10 times those of the actual sample in 
order to reduce the Poisson errors. Using the simulated 
AGNs and the actual flux-area relation of our combined 
sample, we estimated binned SXLFs using three differ- 
ent estimators: £Uf\ N ohs /N mdl (Eq. |) and that of 
PC. The results are compared with the underlying SXLF, 
which was used to generate the simulated AGNs, in Fig. 
[I]. For the models to evaluate A^ mdl , we used the re-fitted 
model using the simulated sample rather than the original 
model. The la errors for the N obs /N lridl (Eq. |) and the 
PC estimators are P oisson errors calculated using Eqs. (7) 
and (12) of Gehrels (|l986|) . On the other hand, the errors 
for the VeT 1 estimator are from Eq. (3) of paper I and 
are inaccurate for bins with a small number of AGNs. 

As shown in Fig. [l| the N ohs /N mdl estimator best rep- 
resents the original model and no estimated point deviates 
from the underlying model by more than 2a. The U a _1 
estimator underestimates the XLF in the lowest luminos- 
ity bin as found in PC. We note that the PC estimator 
also systematically underestimates the LF in this particu- 



lar case of the underlying model and the flux-area relation. 
This is expected because their estimator implicitly builds 
in the assumption — oc const as the underlying LF 
shape. This is much more weighted towards higher lumi- 
nosities than any part of the realistic AGN XLF. Since 
the amount of this bias depends on the underlying model 
and the flux-area relation, as well as the points in the bin 
where the data are plotted, it is not surprising that the 
bias is not apparent in Fig. 2 of PC. They have also com- 
pensated for this bias upon comparing the estimated LF 
with a model. Instead of correcting the estimated binned 
LF using a good model (which our N ohs / N mdl estimator 
does), they calculated the "model-expectated value of the 
estimator" to compare with the estimated value from the 
data. Detailed investigation and comparison of these dif- 
ferent estimators in various cases are beyond the scope of 
this paper. Judging from this simulation, the above dis- 
cussion on biases, and that the exact Poisson errors can be 
used for errors, we choose to use the N ohs /N mdl estimator 
for our plots and tabulation. 



3.4- The binned SXLF results 

Using the N ohs /N mdl estimator, we revised the full SXLF 
plot (Fig. 3 of paper I), as shown in Fig. |[ Instead of 
connecting the data points, we overplotted the analytical 
model for each redshift bin. 

The resulting binned SXLF are listed in Tables ||, |[ 
and [| for different sets of cosmological parameters respec- 
tively. The columns of these tables are — (1) the redshift 
range of the bin; (2) the luminosity range of the bin; (3) 
the number of AGNs in the sample for the bin and the 
number of non-"type 1" AGNs as defined in Appendix A 
of paper I; (4) the number of AGNs expected from the an- 
alytical model derived in Sect. |3.l| . (5) the binned SXLF 
estimated using Eq. || using the model XLF evaluated at 
the central point of the bin, i.e. z — z c (see Table ||) and 
Logi x = (LogL xm i n + LogL xmax )/2, where the subscripts 
min and max signify the borders of the bin in columns (1) 
and (2). The upper and lower errors correspond to Poisson 
errors estimated by Eqs. (7) and (12) of Gehrels ( [1986D 
respectively using S=l (corresponding to the confidence 
of the Gaussian la. When there is no object in the bin, 
the Poissonian 90% confidence upper-limit is given (corre- 
sponding to 2.3 objects). We recommend use of the values 
and errors under this column when, e.g. overplotting ob- 
served SXLF values with model predictions. The ASCII 
versions of these tables with additional columns and sep- 
arate tables with only "type 1" AGNs are provided as a 
part of the source in the preprint archive of this paper 
( |astro-ph/0101279| ). 

Figure || shows that the lowest-redshift, lowest- 
luminosity bin has a significant excess over the two power- 
law analytical expression (AGNs belonging to this bin 
have not been used for the two power-law fit), thus the 
actual underlying SXLF has a much steeper slope than 



Miyaji, Hasinger, Schmidt: ROSAT AGN Luminosity Function 



5 



v 
a 



n in 

o 



10- 4 
10- 5 

10 -6 
10- 7 

10- 8 



^ 10- 9 

>& 

10 -io 

10- 11 

10-18 



I I I I I I I I 1 I I I I I I 

(a) (fi m ,n A ) = (i. 0,0.0) 



+ 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I 



\\\ \\ 



\ 



a z:2.3-4.6 

^ z:1.6-2.3 

• z:0.8-1.6 

O z:0.4-0.8 

■ z:0.2-0.4 

□ z:0.015-0.2 

1 I 1 1 1 1 I 1 1 1 1 



\\\ 4\\ 




_L 



_L 



_L 



(b) (nn A ) = (o. 3,0.0) 



'TV 1^ 



I I I I I I I I I I I I I I I I I I I I I I I I I I I 



42 43 44 45 46 47 
Log L x [h" 2 erg s" 1 ] 



42 43 44 45 46 47 
Log L x [h" 2 erg s" 1 ] 



Fig. 2. The binned SXLFs estimated by Eq.g are plotte d with Poisson errors corresponding to the significance range 
of Gaussian la from Eqs. (7) and (12) by Gehrels ( 1986 ). The data points and error estimates are more accurate than 
Fig. 3 of paper I. Different symbols correspond to different redshift bins as indicated in the lower4eft part of panel 
(a). The symbol attached to a downward arrow indicates the 90% upper limit (corresponding to 2.3 objects) for the 
bin with no AGN in the sample. The best-fit analytical model for each redshift bin in the luminosity range used for 
the fit is overplotted in dashed lines. 



that used for A^ mdl . In order to evaluate the bias caused by 10% lower, thus, the difference is much smaller than the 
this, we have made an N ohs / 7V mdl estimate of this partic- statistical errors for this bin. 
ular bin using the local slope of 7 = 1.7 instead of 7 ~ 0.6 
from the two power-law model. This gave a value about 
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4. Discussion 

We have shown the tables of observed SXLF values for a 
number of standard set of cosmological parameters. These 
values are intended for direct comparison with models and 
plotting with realistic error bars. However, we list a num- 
ber of caveats and sources of uncertainties, and related 
issues. 

— Countrate-to-flux conversion: For the PSPC-based ob- 
servations, where we can limit pulse-height channels, 
the uncertainty in the countrate-flux (in 0.5-2 [keV]) 
conversion is small (± 3% for photon spectral index of 
r = 2.0 ± 0.7. At the faintest end (S x u < 0.5), where 
only the HRI data are available, the conversion rate 
varies by ±40% for the same spectral index range. 

— Optical classification of AGNs: Since different catalogs 
used in this analysis have different criteria for type I 
and type II AGNs, we did not show separate expres- 
sions for these two populations. See appendix A. of 
paper I for the approximate difference in behavior of 
the tentative "type I" sample. 

— Incompleteness: Most of the surveys used in the anal- 
ysis are highly complete or we have selected an appro- 
priate complete subset. In case there is incompleteness, 
we have corrected for it by assuming that the redshift 
distribution/content of the remaining sources are the 
same as the identified ones in the same flux range. This 
assumption is not likely to be the case, considering that 
they have not been identified not because of a random 
cause but because of optical faintness and difficulty in 
obtaining decent optical spectra. The only place that 
this could affect significantly is faintest end of RDS-LH 
(0.17 < Sxii < 0.38), where the identification com- 
pleteness is ~ 80%. (At 5 x i4 > 0.38, the identification 
completeness is <; 95% in any flux range). This could 
affect the behavior of, e.g. the apparent break at the 
low luminosity end in 1.6 < z < 2.3. 

— X-ray spectra/absorption: A serious model composer 
should be aware that the luminosity given here is 
for 0.5 — 2 [keV] in the observer frame. Thus, one 
should compare the model, with their own spec- 



tra^ assumptions (spectral index, absorption, traction ph/0011341) 



— Large-Scale Structure: The lowest redshift bin covers 
0.015 < z < 0.2 and there is some concern about the 
effect of the large-scale structure of the universe, which 
could be confused with the effect of evolution. Zucca 
et al. 1997 found an underdensity of galaxies in the 
local universe out to z~ 0.05. However, it might be 
because of the structure within their survey field of 
~ 27 [deg 2 ] rather than that of the entire space out 
to this redshift. The solid angle surveyed by RBS is ~ 
50% of the sky and the fields of SA-N are scattered in 
various directions. In any event, our lowest redshift bin 
samples a sufficiently large volume of space to z — 0.2 
with a uniform redshift coverage, thus it is unlikely 
that the calculated SXLF is significantly biased by the 
large-scale structure of the universe. 

Acknowledgements. This work is based on a combination of 
extensive R OSA T surveys from a number of groups. Our work 
is indebted to the effort of the ROSAT team and the opti- 
cal followup teams in producing data and the catalogs used in 
the analysis. In particular, we thank K. Mason, A. Schwope, 
G. Zamorani, I. Appenzeller, and I. McHardy for providing 
us with and allowing us to use their data prior to publica- 
tion of the catalogs. TM was supported by a fellowship from 
the Max-Planck-Society during his appointment at MPE. GH 
acknowledges DLR grant FKZ 50 OR 9403 5. We thank the 
referee, T. Shanks, for useful comments. 

References 

Appenzeller I., Thiering I., Zickgraf F.-J. et al. 1998, ApJS 117, 
319 

Avni Y., Bahcall J.N. 1980, ApJ 235, 694 

Bower R.G., Hasinger G., Castander F.J. et al. 1996, MNRAS 
281, 59 

Boyle B.J., Shanks T., Georgantopoulos I., Stewart G.C., Grif- 
fiths R.E. 1994, MNRAS 271, 639 

Comastri A., Setti G., Zamorani G., Hasinger G. 1995, A&A 
296, 1 

Fasano G., Franceschini A. 1987, MNRAS 225, 155 

Fischer J.-U., Hasinger G., Schwope A.D. et al. 1998, Astron. 

Nachr. 319, 347 
Gehrels N. 1986, ApJ, 303, 336 
Gilli R., Risalti G. Salvati M. 1999 A&A 347, 424 
Gilli R., Salvati M., Hasinger G. 2001 A&A in press, (astro- 



of absorbed AGlNs which may depend on lununos- 
ity /redshift), should compute the apparent luminosity 
in the 0.5(1 + z) — 2(1 + z)[keV] range and compare 
it with the values listed in this paper. The latest pop- 
ulation synthesis model based on absorbed and unab- 
sorbed AGNs by Gilli et al. ( |2000| ) has applied this ap- 
proach using the tables shown in this paper. However, 
as discussed in paper I., the no K-correction case corre- 
sponds to a K-correction assuming a power-law photon 
index of T = 2, which is the most representative spec- 
trum for the soft X-ray sources in the sample. Thus 
for many purposes, considering our tabulated values 
as K-corrected SXLF would be accurate enough. 



Hasinger G., 1996 A&AS, 120, C607 
Hasinger G. 1998, AN 319, 37 

Hasinger G., Burg R., Giacconi R., et al., 1993, A&A 275, 1 
Hasinger G., Burg R., Giacconi R., et al. 1998, A&A 329, 482 
Hasinger G., Lehmann I., Giacconi R. et al. 1999 in Highlights 
in X-ray Astronomy in Honor of Joachim Triimper's 65th 
Birthday, MPE Report (Garching:MPE) 
Hasinger G., Altieri B., Arnaud M et al. 2001, A&A, 365, L45 
Hornschemeier A.E., Brandt W.N., Garmire G.P., et al., ApJ 

in press (astro/ph-0042460 
Jones L.R., McHardy I.M., Merrifield MR. et al. 1996, MN- 
RAS, 285, 547 

Lehmann L, Hasinger G., Schmidt M. et al. 2000, A&A 352, 
35 



Miyaji, Hasinger, Schmidt: ROSAT AGN Luminosity Function 



7 



Madau P., Ghisellini G., Fabian A.C. 1994, MNRAS 270, L17 
Marshall H.L., Avni Y., Tananbaum H., Zamorani G. 1983, 
ApJ 269, 35 

Mason K., Carrera F.J., Hasinger G. et al. 2000, MNRAS 311, 
456 

McHardy I., Jones L.R., Merrifield M.R. et al. 1998 MNRAS 
295, 641 

Miyaji T., Ishisaki Y., Ogasaka Y., Ueda Y. et al. 1998, A&A 
334, L13 

Miyaji T., Hasinger G., Schmidt M. 1999, in Highlights in 
X-ray Astronomy in Honor of Joachim Trumper' s 65th 



Table 3. The full binned SXLF values - 1. 



Bir thday, M PE Report (Garching:MPE) in press, (astro 
ph/|9809398|, M99a) 



Miyaji T., Hasinger G., Schmidt M. 2000, A&A 353, 25 (paper 
I) 

Miyaji T., Hasinger G., Schmidt M. 2000, Adv. Sp. Res. 25, 
827 

Mushotzky R.F., Cowie L.L., Barger A. J., Arnaud K.A. 2000, 

Nature 404, 459 
Page M.J., Mason K.O., McHardy I.M., Jones L.R., Carrera 

F. J. 1997, MNRAS 291, 324 
Page M.J., Carrera F.J. 1999, MNRAS 311, 433 
Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P. 

1992 Numerical Recipes in Fortran (Cambrideg: Cambridge 

Univ. Press), 640 
Schmidt M. 1968, ApJ 151, 393 
Schmidt M, Green R.F. 1983, ApJ 269, 352 

the Universe, MPE Report 263 (Garching:MPE), 395 
Schmidt M., Hasinger G., Gunn J. et al. 1998, A&A 329, 495 
Schwope A.D., Hasinger G., Lehmann I. et al. 2000, Astr. 

Nachr. 321, 1 

Voges W. 1994, in Basic Space Science, Haubold H.J. & Onuora 
L.I. (eds.) 202 

Voges W., Aschenbach B., Boiler Th. et al. 1999, A&A 349, 
389 

Wisotzki L. 1998, Astr. Nachr. 319, 257 

Zamorani G., Mignoli M, Hasinger G. et al. 1999, A&AS 346, 
732 

Zickgraf F.-J., Thiering I., Krautter J. et al. 1997, A&AS, 123 
103 

Zucca E., Zamorani, G., Vettolani, G. et al. 1997, A&A 326, 
477 



LogL x a 



,0a) = (1 
jV obs 



0,0.0) 

jymdl 



(1) 


(2) 


(3) 


(4) 


.015-0.2 


41.30-41.70 


5(3) 


1.0 


.015-0.2 


41.70-42.40 


8(2) 


6.3 


.015-0.2 


42.40-43.00 


23(3) 


26.8 


.015-0.2 


43.00-43.50 


61(7) 


58.4 


.015-0.2 


43.50-43.80 


56(3) 


52.6 


.015-0.2 


43.80-44.20 


55(4) 


65.5 


.015-0.2 


44.20-44.50 


35(1) 


31.6 


.015-0.2 


44.50-45.00 


30 


25.0 


.015-0.2 


45.00-45.70 


1 


3.0 


0.2-0.4 


41.70-42.10 


1 


0.8 


0.2-0.4 


42.10-42.70 


3(2) 


4.3 


0.2-0.4 


42.70-43.30 


18(5) 


15.1 


0.2-0.4 


43.30-43.80 


27(3) 


26.3 


0.2-0.4 


43.80-44.30 


25(3) 


29.0 


0.2-0.4 


44.30-44.90 




23.5 


0.2-0.4 


44.90-45.40 


13 


11.3 


0.4-0.8 


42.30-42.70 





1.1 


0.4-0.8 


42.70-43.30 


13(7) 


10.8 


0.4-0.8 


43.30-43.60 


14(4) 


11.7 


0.4-0.8 


43.60-44.20 


38 


43.2 


0.4-0.8 


44.20-44.80 


19(1) 


16.8 


0.4-0.8 


44.80-45.40 


10 


8.9 


0.4-0.8 


45.40-46.50 


5 


6.5 


0.8-1.6 


42.70-43.30 


5 


2.9 


0.8-1.6 


43.30-43.90 


23 (2) 


26.9 


0.8-1.6 


43.90-44.50 


55 


53.0 


0.8-1.6 


44.50-45.10 


39 


36.0 


0.8-1.6 


45.10-45.70 


5 


9.0 


0.8-1.6 


45.70-46.20 


4 


3.1 


0.8-1.6 


46.20-46.90 


4 


2.3 


1.6-2.3 


43 60-44 20 


g 


13.1 


1.6-2.3 


44.20-44.80 


18 


12.9 


1.6-2.3 


44.80-45.50 


14 


10.5 


1.6-2.3 


45.50-46.10 


2 


2.3 


1.6-2.3 


46.10-46.80 





1.5 


1.6-2.3 


46.80-47.40 


1 


1.0 


2.3-4.6 


43.70-44.10 


2 


3.4 


2.3-4.6 


44.10-44.80 


12 


10.5 


2.3-4.6 


44.80-45.40 


7 


4.3 


2.3-4.6 


45.40-46.20 


2 


3.2 


2.3-4.6 


46.20-47.00 





1.4 


2.3-4.6 


47.00-47.50 


2 


1.0 



d Log _L X 

(5) 



231^ 
2.8±i;g) 

7.7ii:») 
3.9±8:i) 

+0.2\ 
0.2) 
+0.6N 
0.6) 

1.3 4 " 



1.6 
4.1 

+0.2 
0.2 
2 1+ 0.4 

3.3; 



0.4/ 
6.8 N 
2.7) 



1.1 

2.8 
1.9 
5.9 
1.1 
9.6 
4.7 



< 7.9 
3 

°- ± -0.9 
+0.6 
0.5 

3.9+ 07 



1.81 



3.3 
1.5 
1.3 



8.9Z 



+5.6 

-3.8 

2.4i°; 6 5 

5.6+ 1 ' 



1.1 
1.2 
6.1 



1.82 



+2.3\ 
0.9^ 
+2.5\ 
1.5) 
+0.5\ 
OA) 
+1.3) 

+0.3-, 
0.2) 
+2.2-, 
1.9J 
+ 1.6-, 
1.3) 



■ 10" 


-4 


■ 10" 


5 


■ 10" 


-5 


■ 10" 


-6 


■ 10" 


-6 


■ 10" 


-8 


■ 10" 


-9 



0.6; 

+0.9> 
0.7) 
+0.6> 
-0.5) 
0.8\ 
0.6) 



0.9) 
+0.7 \ 
0.5) 
+0.9-, 
-0.6) 
4.5-, 
2.9/ 



+0.8- 

0.6, 

5.8±i;|: 

9 c+0.8< 

6.7±S;g: 
< 4.6- 



10" 


io 


■ 10" 


-5 


■ 10" 


-5 


■ 10" 


-6 


• 10" 


-7 


■ 10" 


-8 


■ 10" 


-9 


10" 


11 


• 10" 


-5 


■ 10" 


-6 


■ 10" 


-7 


■ lQ- 


-9 



10" 



(3 



+4.3\ 
2.3) 
+2.2-, 
1.7) 
+2.3\ 
1.7) 
+9.4-, 
5.0/ 

< 5.4- 10~ 11 

tr 5 ) ■ 10- 



3.6 
6.0 
4.6 
7.8 



Notes: a L x [ftJ 2 ergs" 1 ] in 0.5-2 [keV]. b [hl Mpc" 3 ]. c The 
AGNs in this row are outside of the luminosity range used 
for the model fit. 



8 Miyaji, Hasinger, Schmidt: ROSAT AGN Luminosity Function 

Table 4. The full binned SXLF values - 2. Table 5. The full binned SXLF values - 3. 
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(3.8+?-J) 


■ 10" 


-7 
















2.3-4.6 


45.40-46.00 


7 


4.6 


(4.2t? : £) 


■ 10" 


-8 


2.3-4.6 


44.10-44.80 


7 


9.0 


(4.6ir 7 ) 


■ 10" 


-6 


2.3-4.6 


46.00-46.60 


1 


1.9 


(i-ol2:S) 


■ 10" 


-9 


2.3-4.6 


44.80-45.40 


8 


6.6 


(4.H1:2) 


■ 10" 


-7 


2.3-4.6 


46.60-47.10 





1.2 


< 3.2 ■ 10" 10 


2.3-4.6 


45.40-46.20 


7 


5.8 


(i-9±8:?) 


■ 10" 


-8 


2.3-4.6 


47.10-47.60 


2 


1.0 


(3.71^) 


10" 


11 


2.3-4.6 
2.3-4.6 


46.20-47.00 
47.00-47.63 


1 

2 


2.2 
1.3 


(2.iil: 7 ) 
(3.ol 3 :^) 


• 10" 

• 10" 


10 
11 



7-8l 6 5 ) 



1.2" 1 



1.4x 
0.8 / 
1 o+0.4n 
± -°-0.3/ 

4 , +i.o 



1.1 



0.9; 

+0.2s. 
0.2) 

1- 3±8:l) 

2- Oll'f,) 



• 10"° 

■ 10" 5 

■ 10" 5 

■ 10" 6 

■ 10" 6 
■10" 7 

■ 10" 8 

• 10 



10 



Notes: a L x [/i^" 2 erg s" 1 ] in 0.5-2 [keV]. b [/i| Mpc" 3 ]. c The 
AGNs in this row are outside of the luminosity range used 
for the model fit. 



Notes: a L x [/ig 2 ergs" 1 ] in 0.5-2 [keV]. b [hf >0 Mpc" 3 ]. c The 
AGNs in this row are outside of the luminosity range used 
for the model fit. 



